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Abstract 

Solving the QCD renormalization group equation at the 2-loop and 3-loop orders we obtain 
exphcit expressions for the coupUng as a function of the scale in terms of the Lambert 
W function. We study the nature of the "Landau singularities" in the complex plane 
and show that perturbative freezing can lead, in certain cases, to an analyticity structure 
that is consistent with causality. We analyze the Analytic Perturbation Theory (APT) 
approach which is intended to remove the "Landau singularities", and show that at 2- 
loops it is uniquely defined in terms of the Lambert W function, and that, depending on 
the value of the first two /3 function coefficients /3o and it is either consistent with 
perturbative freezing (for /3i < — /Jq) with an infrared limit of — /3o//?i or leads to a non- 
perturbative infrared coupling with a limit of l//3o (for /?i > — /?o)- The possibility of 
a causal perturbative coupling is in accordance with the idea that a purely perturbative 
Banks-Zaks phase with an infrared fixed-point exists in QCD if the number of flavours 
(A^/) is increased. The causality condition implies that the perturbative phase is realized 
for > 10. 
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1 Introduction 



Due to asymptotic freedom [1], physical quantities in QCD at large momentum 
transfers = — g^, where is a space-like momentum-squared, can be calculated 
as power expansions in the coupling constant ^(Q^) = as{Q'^)/'K. This seems as a 
reasonable expansion^ since the coupling vanishes at this limit according to: 

where 

/i^o = ^(ll-^iV/) (2) 

and A is the QCD scale. x{Q'^) in (1) is a solution to the 1-loop renormalization 
group (RG) equation 

Going to lower Q^, x{Q'^) becomes larger, and higher loops have to be taken into 
account in (3), as well as in the expansions that describe physical quantities in terms 
of x{Q'^). The RG equation at the (n -|- l)-th loop order is 

dx 

<5'l7^ = l^i.^) = -Po + CX + C2X^ + ■■■ CnX""). (4) 

For Q"^ below A^, non-perturbative effects, which are non-expandable in x{Q'^), 
become dominant, and the perturbative expansion becomes useless. However, in the 
intermediate regime, the perturbative solution can be fitted to the data, provided 
it is supplemented by power-like terms. These terms are non-perturbative but they 
can be characterized by perturbation theory, as they are related to the large order 
asymptotics of the perturbation series, and in particular to renormalons (see [2] and 
refs. therein). It is the non-convergence of the perturbative expansion, that provides 
some information on the non-perturbative terms. 

Another indication that the perturbative result is incomplete, and cannot de- 
scribe the low physics unless it is supplemented by non-perturbative corrections, 
comes from considering its analyticity structure in the complex plane: a generic 
QCD observable, that depends on a space- like momentum variable Q^, is expected 
to be an analytic function of in the entire complex plane, except the negative 
real (time-like) axis. Singularities on the time-like axis are meaningful since they 
correspond to production of on-shell particles, while existence of singularities away 
from the time-like axis would violate causality. Thus causality constrains the func- 
tional dependence of observables on Q^. The analyticity condition is equivalent to 
the requirement that x{Q^) obeys the following dispersion relation: 

Puis) = --lm{x{-s)} (5) 

TT 



^Eyen though the series actually does not converge, and in-fact it is non Borel-summable. 
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with 



x{Q') = - / 

JO 



oo 



ds 



(6) 



where I3ji{s)^ with < s < oo, is called the spectral density function. 

A 1-loop perturbative result for a generic QCD observable depends on the cou- 
phng (1) and therefore contains a "Landau-pole" at — Is?. This pole is located 
on the positive real axis (the space-like axis) and therefore it is non-physical. This 
pole is expected to be washed out when further perturbative and non-perturbative 
corrections are added. In general, higher loop perturbative results for the coupling 
have a more complicated analyticity structure which is still inconsistent with causal- 
ity. Thus, in general, causality can by saved only be inclusion of non-perturbative 
terms. 

A special case in which the perturbative result by itself can be consistent with 
causality, is when the perturbative j3 function (4) has a zero. Then, the coupling 
reaches a finite value in the infrared limit and thus it is finite for any real and 
positive Q^. To be consistent with causality, the coupling should be non-singular 
in the entire plane except the time-like axis, and therefore there may be cases 
where freezing occurs but causality is still violated, due to complex singularities. 

In real-world QCD, with only three light flavors, it seems from the flrst few 
terms in the j3 function series that there is no perturbative freezing [3]. In this 
case, the standard perturbative approach always faces the "Landau singularity" 
problem. Refs. [4, 5] revive a natural solution to this problem through the so- 
called Analytic Perturbation Theory (APT) approach. The APT coupling was used 
in [7] for discussing power corrections within the more general (non-perturbative) 
dispersive approach of [6], and was also considered there (as well as in refs. [2, 8]) 
as a possible model for non-OPE power terms. 

In this paper we solve the 2-loop and 3-loop RG equation giving explicit expres- 
sions for the coupling as a function of the scale in terms of the Lambert W function. 
This enables us to study the location and the nature of the "Landau singularities" 
in the complex plane. In particular we flnd exact criteria which determine when 
perturbative freezing a la Banks-Zaks [10] leads to an analyticity structure that 
is consistent with causality. Using the Lambert W solution, we then analyze the 
APT approach and show that it is either consistent with perturbative freezing (for 
c < — /^o < 0) with an infrared limit of — 1/c, or leads to a non-perturbative coupling 
with an infrared limit of l//3o (for c > — /5o). 



2 The Lambert W function — exact explicit cou- 
pling at 2-loops 

We start with a 2-loop RG equation 



(3{x) 



dx 



Pox^ (l + cx) 



(7) 



2 



where t — ln((5^/A^). Note that the value of A here should be different from the 
one in (1), so that the same phenomenological value for will be obtained. 

When higher-order corrections are added, A should be further adjusted, but for 
simplicity we use the same notation throughout. The 2-loop coefficient c, which is 
renormalization scheme independent, is given by: 



102 - -Nf 
3 \ 



(8) 



A straightforward integration of (7) for the asymptotically-free case {Po > 0) yields: 

(9) 



/?oln(gVA2) = i-cln 

X 



1 

- + c 

X 



At this stage it is already clear that x{Q'^) has a cut on the negative real axis, being 
a function of ln(Q^/A^). The problem of inverting (9) makes it difficult to study the 
singularity structure. In [4], for instance, the inversion of (9) relies on the assumption 
that the 1/x term is much larger than the logarithmic term. Although correct in 
the deep perturbativc region, this approximation is inadequate for studying the 
singularity structure of x{Q'^). In particular, this approximation does not allow 
perturbative freezing at all, since in the latter case the logarithmic term becomes 
dominant in the infrared region. 

Luckily, an explicit solution of the 2-loop RG equation (7), i.e. an inversion of 
(9), can be written in terms of the so-called Lambert W function [11]. W{z) is 
defined by 

W{z)eycp\W{z)\^ z (10) 

The coupling is then given by: 



c l + W{z) 




(11) 

-/3o/c 



W{z) is a multi-valued function with an infinite number of branches, denoted by 
Wn{z). We follow [11] as for the division of the branches and notation. The require- 
ment that x{Q'^) is real and positive for a real positive Q"^ (at least for Q"^ » A^), is 
sufficient to determine the physical branch depending on the sign of c: 

(a) for c > 0, z < and the physical branch is W_i{z). This branch is a real 
monotonically decreasing function for z e (— 1/e, 0), with W-i{z) e (— oo, —1) 
(outside this range W-i{z) is complex). Thus, the ultraviolet limit corresponds 
to z — 0~, W-i{z) — — cxD, and x — O"*", as required by asymptotic 
freedom. In the infrared region, below the "Landau singularity" (see further 
discussion later) x is complex. 
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(b) for c < 0, 2; > and the physical branch is the principal branch, Wo{z). This 
branch is a real monotonically increasing function for z G (— l/e,oo), with 
Wq G (— 1,00) (outside this range, Wo{z) is complex). For -2 > 0, Wo{z) G 
[0,00). The ultraviolet limit corresponds to z — > 00, Wq{z) — > 00 and 
X — 0"^. The infrared limit corresponds to z — 0"*", Wo{z) — > 0"*", and 
X — > (— l/c)~ which is consistent with the Banks-Zaks perturbative fixed- 
point value xpp — — 1/c. 

The solution in (11) can be quite useful whenever a two- loop evolution of the 
QCD coupling is required. It clearly yields more accurate results than the standard 
expansion of the coupling in 1/ ln{Q'^/ A^) . Note that the latter can be obtained from 
(11) using the asymptotic formulae for W{z): for c < one uses the asymptotic 
expansion of Wo{z) at large positive z that starts with Wo{z) ~ In(^) — ln(ln(2;)), 
and for c > one uses the asymptotic expansion of W-i{z) at small negative z that 
starts with W^i{z) ~ ln(— 2;) — ln(— ln(— 2;)). 

Other possible applications of (11) arc in the context of resummation methods, 
such as the iterated construction suggested by Maxwell in [12] and the improved 
Baker-Gammel approximants suggested by Cvetic in [13] which both generalize the 
diagonal Pade approximants approach [14]. 

The general idea in these resummation methods is that knowledge of the first 
few orders in a perturbative expansion, together with the RG equation, can be used 
to construct approximants to physical quantities that have a better accuracy than 
the truncated perturbative series, and exhibit reduced renormalization scale and 
scheme dependence (for related ideas, see [15, 16]). Diagonal Pade Approximants, 
as opposed to truncated series, are exactly invariant with respect to an arbitrary 
renormalization scale transformation, so long as the transformation respects the 1- 
loop evolution form (i.e. an evolution related to the 1-loop f3 function) [14]. In 
order to achieve scale invariance beyond this approximation, it turns out [12, 13] 
that one has to use more complicated functions, which are related to the inverted 
solution of the some higher-order RG equation. Using the explicit inverted solution 
of the 2-loop RG equation (11) the methods of [12, 13] can be easily implemented 
in practice and also their mathematical structure can be further studied. 

It is worth mentioning that another way-^- to invert (9) was suggested in the past 
by Khuri and Ren [17]: they wrote the 2-loop couphng in terms of -F(C), which is 
defined by 2F(C) — exp[F(C)] -|- 1 = C- The latter equation was considered earher 
in different physical contexts (sec rcfs. in [17]). There is, of course, a one-to-one 
correspondence between the Lambert W solution of (11) and that of [17]: 



--exp[F(C)] = w(--exp 
In the following we stick with the Lambert W solution. 



i(C-l) 



12) 



^An exact solution of (9) in term of a (modified) Borel representation can also be found in [18]. 
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3 The singularity structure of the 2-loop coupHng 



Having found an explicit solution for coupling (11), we would like now to analyze its 
singularity structure in the complex plane and, in particular, to find when it is 
consistent with causality. To do so it is essential to define the analytical continuation 
of x{Q'^) to the whole complex plane, namely to specify the branch of W{z) that 
is used in (11) for a generic = \Q^\e'^'^, where —n < (j) < n. 

At this stage it is necessary to give a brief description of the singularity structure 
of the Lambert W function; more details can be found in [11]. The partition of the 
complex W plane between the different branches, Wq, W±i and W±2., is shown by 
dashed lines in fig. 1 (the other lines in fig. 1 will be discussed later). It is important 
to realize that the singularity structure of the different branches Wn{z) is different. 
Wo, W^i and Wi share a common branch point at 2; = — 1/e, at which W = —1. 
The z = — 1/e cut in all three branches is chosen to be 2; G {—00,—l/e). While 
for Wo, z — —1/e is the only singularity, other branches Wn{z) with n 7^ have a 
branch point at 2; = with a cut at (—00, 0). The z ~ cut is the only singularity 
in Wn for n > 2. Note that W{y) obeys the following symmetry [11], 



and therefore it is possible to obtain W{z) and x{Q^) on the upper half-plane from 
those on the lower half-plane. Note also that the branches W_i{z) and H^i(^) have a 
common real limit along the cut for z e (—1/e, 0). It is only a matter of convention 
that the negative real axis (W e (—00, —1/e)) belongs to rather than to W^. 

The criterion we use for defining the analytical continuation is that for \Q'^\ ^ A^, 
the coupling x(Q^), and therefore also W{z), will be continuous as a function of the 
phase of Q^. From the analyticity structure of W{z) described above it is clear that 
as long as the phase of z does not reach ±7r for any Q"^ in the first sheet, the only 
relevant branch is the physical one, i.e. the one that yields a real positive x for real 
positive Q"^: Wq for c < and W±i for c > 0. Let us define z = \z\e^^. Using (11) 
the condition — tt < 5 < tt can be directly translated to conditions on /?o and c, and 
therefore further division of the c axis is suggested as follows: 

(a) for c < —Po < we have S = — (/9o/c)0, and for G {—71,71) we always have 
— TT < — |/9o/c|vr < S < |/9o/c|7r < tt. Thus Wo{z) is the only relevant branch, 
its boundary is never reached and no singularity is encountered. We conclude 
that here the coupling has a perturbative fixed-point at xpp = — 1/c, and an 
analyticity structure that is consistent with causality. 

(b) for —(3o < c < 0, we find that 6 = — (/3o/c)0 reaches ±7r at ±0o, with 



Thus, when has a phase of ±0o, the boundary of the Wo{z) branch is 
reached. Consequently, an image of the cut z G (—00, —1/e) appears on the 



Wl^iy*) = W^{y) 



(13) 



00 = \c/(3q\ tt. 



(14) 
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first sheet in the Q"^ plane. The branch point corresponding to z — — 1/e 
appears at 

where 

g2 = a2|c|--//3o. (16) 

The cuts in the plane and the analytical continuation for |0| > 0o will be 
discussed later. 

(c) for c > /9o/2, the solution is given by W-i{z) with 5 = +7r — {Po/c)4> for > 
and by Wi{z) with 5 = -tt- (/5o/c)0 for < 0. For G (0,7r), 6 > {l-po/c)7i 
and therefore it never reaches the W^i branch boundary at 5 = — tt. Similarly, 
for (p G (— TT, 0), 6 < (— l + /3o/c)7r and therefore it never reaches the Wi branch 
boundary a,t 5 — +7r. Thus, the only relevant branches are W±i and the only 
singularity encountered is the one at 2; = —1/e with the cut z G (—00, —1/e) 
that separates between Wi and W^i. From (11) wc find that the z = —1/e 
singularity appears on the positive real axis, at = Q^, where Qq is given 
in (16). The cut z G {—00, —1/e) corresponds to a cut on the infrared section 
of the positive real axis, G (0, Ql). 

(d) for < c < Po/2, 6 = ±7r — (/3o/c) 4> and the boundary of W_i at 5 = —tt, 
and that of Wi at S = tt, is reached at = ±0i, where 0i = 2 (c//?o) tt. Like 
in case (b) above, it is required to define W{z) for > 0i. This will be done 
soon. 

In the two cases with large \Po/c\, (b) and (d) above, we found that the boundary 
of the physical branch is reached, and therefore a definition of W{z) in eq. (11) 
for 101 > 00, in (b), and |0| > 0i, in (d), is required. The criterion that W{z) 
should be continuous as a function of for large \Q'^\ implies a particular definition 
for the analytical continuation of W{z): starting with a given branch iy|„|(2;) and 
increasing |0|, then when the branch boundary is reached one switches to the next 
branch l^|„|+i(z). 

Let us illustrate the above starting with the case < c < /?o/2, (d) above, and 
consider a gradual increase in |0|, the phase of Q^, in the lower half-plane < 0. 
For a positive real Q^, z < and W{z) is just on the boundary between W^i 
and Wi, shown by the dashed line in fig. 1. For a small negative 0, z is slightly 
below the negative real axis {S — —tt — (/3o/c)0) and the solution is within the 
branch Wi{z). This is exemplified in fig. 1 by the case (7 = 0.1 (0 = — O.Itt). 
When = —01 = — 2(c//5o)7r, 6 = tt, z becomes negative real again, and so the 
cut z G (—00,0) which is the upper boundary of the Wi{z) branch is reached. 
Continuity of W{z) is obtained if, and only if W2{z) is used for < — 0i. In fig. 1, 
this is exemphfied by the case g — 0.6 (0 = — 0.67r). Note that in the specific 
example given in fig. 1, the time- like axis = —tt is reached within the W2 branch 
(the dot-dash line), but in general, depending on the ratio Po/c, higher branches of 
W may become relevant. 
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Next, consider the case — /3o < c < 0, (b) above, where the positive real 
solution is the positive real W axis in the Wq branch, described by the g = line 
in fig. 2. For a small negative (p, 6 = — (/3o/c)0 is still far from the 6 = —n limit 
and the solution is given by Wq{z). In fig. 2, this is exemplified by the case g = 0.1 
{4> = — O.Itt). For (f) = —00 = (c//3o)7r, the cut ai z — (— oo, — 1/e) which is the 
lower boundary of the Wq branch is reached. Continuity of W{z) at large (the 
right side of fig. 2) implies that for (f) < — the branch W^i{z) should be used. 
In fig. 2 this is exemplified by the case g = 0.7 (0 = —O.Tn). The possibility that 
higher branches {W\n\ for n > 1) will become relevant (depending on the ratio Po/c) 
exists also here. Just hke in the previous case, these branches are reached through 
the z — cut. 

A unique feature of the analytical continuation we perform appears in case (b) 
(see fig. 2): this is a "phase transition" from the Banks-Zaks non-trivial infrared 
fixed-point to the trivial one, that occurs at = ±0o- For any \4>\ < 4>q (the lines 
with g < 0.6 in fig. 2) W{z) flows to zero in the infrared, which corresponds to 
the perturbative flxed-point at xpp = — 1/c. On the other hand, for any |0| > 0o 
(the lines with g > 0.7 in fig. 2) |Vr(2;)| becomes infinite in the infrared, which 
corresponds to x^Q"^) — > 0. The boundary between the Wq branch and the W^i 
branch in fig. 2 separates between the two regimes. A fiow from the ultraviolet 
down to the infrared for (j) = (j)Q {4> = —{2/3)n in fig. 2) leads to a fork at 1^ = — 1 
(corresponding to — Qq) from which there are two ways to continue towards 

= 0, one is to the right, within the Wq branch leading to W = 0, and the 
other to the left, on the boundary between the W^i and Wi branches, leading to 
W = —oo. In other words, choosing the analytical continuation beyond |0| = 0o 
as we did, guarantees continuity for G {Ql, oo), but leaves a discontinuity in the 
complex plane, along the = ±0o directions for any e (0, Ql). 

Let us now summarize how we choose the branch of the Lambert W function in 
(11) according to the proposed analytical continuation of the coupling to the entire 
complex plane. It is convenient to determine the branch from d = — (/3o/c) 0, 
where is the phase of Q^, as before^. Given d, we use the branch Wn such that for 
c < 0, 

de{{2n-l)n, {2n+l)n] (17) 

and for c > 0, 

d e {2{n - 1) TT , 2n7r] n>l (18) 
(i e (2n7r, 2(n+ l)7r] n < -1 

For instance, for c < 0, W^i is used for d e (— Stt, — tt] and Wq is used for d e (— tt, tt], 
as shown in fig. 2. For c > 0, VFi is used for d e (0, 27r] and W2 is used for 
d G (27r, 47r], as shown in fig. 1. 

Note that the lower half of the complex plane (0 < 0) always corresponds to a 
positive imaginary part for the coupling. For c < the lower half-plane corresponds 

^Note that d is not the phase of z, since the latter has an additional n term for c > 0. 
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to d < which, according to (17), refers to n < i.e. Im{l^} < that yields 
lm{x{Q^)} > in (11). For c > 0, < corresponds to d > which according to 
(18) refers to n > 1 i.e. Im{W^} > that again yields lra{x{Q^)} > in (11). 

We are now in a position to address the question we started with, namely what is 
the singularity structure of x{Q'^) in the complex plane. There are three different 
possibilities, depending on the values of c and /3o > 0, which are shown in fig. 3^: 

(a) for c < —Po < the 2-loop coupling has a non-trivial perturbative infrared 
fixed point at xpp — —1/c and an analyticity structure that is consistent with 
causality. Note there is no pole in the denominator of (11), since whenever 
W is real, it is positive. From considering only the 2-loop (3 function it is not 
possible to exclude the physical relevance of this fixed-point, and there is no 
indication of the existence of non-perturbative effects. 

(b) for — /3o < c < the j3 function seems to have a fixed-point, as in a), but 
since xpp = —1/c is large here, this fixed-point is probably not reliable - it 
can be washed out by higher-loop contributions or by non-perturbative effects. 
Existence of the latter is strongly suggested by the presence of causality vio- 
lating singularities: there are two branch points at Q"^ ~ Q'^e^'^'^° with radially 
oriented cuts that end at = 0. The pole in the denominator of (11) at 
W — —1 falls right on top of these branch points. 

(c) for c > there is no infrared fixed-point. The singularity structure violates 
causality, due to the branch point on the space-like axis at = Qo- The pole 
in the denominator of (11) at = —1 falls right on top of this branch point. 
Note that for c > it is not important whether /3o/c is large or small: in any 
case the only singularity is the single image of the z — — 1/e branch point, 
the starting point of the cut that separates between Wi and W-i (There are 
no similar singularities for higher branches \n\ > 2 of the W function). The 
non-physical cut on the space-like axis is probably removed from any physical 
quantity by non-perturbative effects. 

The truncated 2-loop (3 function (7) can be considered as a legitimate choice of 
renormalization scheme - the so-called 't Hooft scheme. However, in certain cases 
this may be a peculiar choice of scheme for calculating observable quantities: the 
truncation of the (3 function at the 2-loop order seems quite arbitrary. Therefore, 
it is interesting to see to what extent our solution depends on this specific choice of 
scheme. 

^Some of these conclusions have been anticipated by Uraltsev [19]. 
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4 The Lambert W solution at 3-loops 



Given the above motivation we would like to generalize our results by considering a 
generic 3-loop RG equation: 

I3(x) = ^ = -Pox'' {1 + CX + C2x2) (19) 

where again t = ln((5^/A^). A straightforward integration of (19) yields a function 
that involves arctanh terms in addition to the terms of the form appearing in (9). 
Therefore, it is even harder to explicitly invert this function. To avoid this difficulty, 
we use the following trick: we alter the 3-loop P function slightly, by taking its x^[l/l] 
Pade Approximant (PA)II: 



/?P.(x) = (20) 

1 - [Co C)X 



This change does not limit generality as long as we restrict our interest to 3-loops, 
since the difference between the /5 function in (20) and that in (19) is of the order 
C(x^). Indeed, using (20) instead of (19) is expected to change the singularity 
structure of x{Q'^), but we insist that (20) corresponds to a choice of scheme at 
4- loops and beyond which is just as arbitrary as the truncation in (19). Using (20) 
we obtain: 

/?oln(gVA2) = --cln^^ ■ - 



X 



+ c 

X 



(21) 



and finally. 



c l-{c2/c^) + W{z) 

(22) 

1 



— exp 



-1 + (C2/C') - Pot/C 



Just like in the 2-loop case (11), the sign of z and therefore also the physically 
relevant branch of W{z) are determined according the sign of c: for c > 0, z < 
and the physical branch is W_i{z), taking real values in the range {—oo, —1), while 
for c < 0, 2; > and the physical branch is Wo{z), taking real values in the range 
(0,00). 

Note, that the only significant difference between the 3-loop solution (22) and 
the 2-loop solution (11) is in the relation between W{z) and x{Q^). This is because 
the difference in the definition of z can be swallowed in a redefinition of the scale 
parameter 

A2 ^ A2 ^ A2eC2/(/3oc)_ (23) 

Consequently, we immediately know almost everything concerning the singularity 
structure of the PA-improvcd 3-loop coupling: it has the same type of branch points 
and cuts that are described in fig. 3. The difference is, however, that the simple 

IIPA's were found useful in various applications in perturbative QCD - see refs. [14, 20]. 
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pole due to a zero in the denominator of x{Q'^) in (22) at W{z) — 02/0^ — 1, is not 
obtained on top of the branch point (s) &i W — —1 (which corresponds instead to 
the pole of the PA improved j3 function). Let us consider the three possible cases: 

(a) If c > then when W is real {W — VF_i) it obeys W{z) < — 1 and then any 
C2 < would result in a pole on the space-like axis. 

(b) If c < -(5q then when W is real {W = Wq) it obeys W{z) > and then if 
C2 > the pole will appear on the space-like axis (for C2 < (? the coupling is 
causal) . 

(c) If — /3o < c < 0, 1^ obtains (or approaches) any real value: W — Wq G (0, 00) 
on the space-like axis, W — > G (-1, 0) for |0| — > 0o and W — > W±i e 
(—00, —1) for 101 — ^ 00 . Therefore for —[3q < c < 0, the pole will always 
appear, either on the space-like axis (if C2 > c^), or at |0| = 0o (if C2 < c^)- 

We conclude that the singularity structure of the PA-improved 3-loop coupling 
is not much different from the 2-loop coupling: there are regions of the parameter 
space for which perturbative freezing is consistent with causality. In particular, this 
is true for small enough and c < 0, which is the starting point for the Banks-Zaks 
expansion. In general, however, higher-order terms in the /3 function create a more 
complicated singularity structure in the infrared region, which is inconsistent with 
the analyticity requirement. 

5 The Analytic Perturbation Theory approach 

Recently it was suggested [4] to construct a causal coupling constant by analytically 
continuing the coupling to the time-like axis (i.e. looking at x{Q^ — —s) with s > 0), 
use the discontinuity, or equivalently the imaginary part, to define a spectral function 
I3r{s) as in (5) and finally use the dispersive integral (6) to construct a new space-like 
coupling 2;apt(Q^)- APT stands for "Analytic Perturbation Theory" since in this 
approach x^p^{Q'^) has the required analyticity structure by construction, being 
defined through the dispersion relation (6). 

The authors of refs. [4, 5] show that the APT couphng has some remarkable 
features such as a universal infrared limit of l/,/?o which is almost independent of 
the renormalization scheme and of the order of the f3 function one starts with. A 
general argument for the universality of the APT infrared limit (for c > 0) was also 
given in [7]. With the explicit expression for the 2-loop coupling in (7), we can now 
directly check these conclusions. 

Before dealing with the 2-loop case, let us briefiy review the APT results at 
1-loop, where both the spectral function and a^^pT can be easily obtained in a 
closed form. One starts with the 1-loop (3 function (1) and performs an analytical 
continuation by substituting Q'^ — —s — ie where s > and e > is small. In the 
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limit e — > one obtains, 

Poll 



lm{x{—s)} 



(/3o7r)2 + (Pot) 



2 



where t = ln(s/A^). The spectral density is 

Pn{t) = -hm{x{-s)} = - ] (25) 

TT /5o [(ln(s/A2))^ + 7r2j 

Finally, one constructs the corresponding space-like effective coupling, through the 
dispersion relation (6). The integral can be performed analytically, and yields 



1 

+ 



ln(Q7A2) A2-g2 



(26) 



The first term in (26) is just the 1-loop perturbative result and the second term 
exactly cancels the "Landau-pole". Since the second term is a power correction, 
it does not alter the ultraviolet behavior. By construction, xji^p^{Q'^) has a cut 
at Rc{Q2j < and no other singularities in the complex plane, and is therefore 
consistent with causality. The coupling has a non-perturbative infrared fixed-point 
at xapt(O) = 1/ Po- 
ll is clearly of interest to see how these APT results change at 2-loops. Since we 
already know how to analytically continue the space-like coupling to the entire 
plane, we can now examine x^Q"^) on the time-like axis. Before doing so, we stress 
once more that for the uniqueness of the APT solution it is necessary to impose 
a continuity condition on x{Q'^) for large \Q'^\. This condition is implemented by 
starting from x{Q^) on the space-like axis with ^ A2, and demanding continuity 
of a:((5^) as the phase of is changed. This prescription is essential, since the naive 
alternative of "taking a shortcut" and going directly to the time-like axis by trying 
to invert (9) for negative leads to an ambiguity in the choice of the branch of W 
and in the corresponding APT spectral function. 

Actually, perturbatively speaking, there can be no ambiguity in the definition of 
x(Q^) at complex Q"^. To see this one first writes the solution of cq. (4) in the stan- 
dard "non-improved" way as x^Q"^ / fi^,Xo), where xq = x{Q'^ = fi"^), and expands in 
powers of xq. The resulting coefficients are polynomials in log((5^/ A*^) (the standard 
RG logs) and can all be expressed in terms of the P function coefficients. If the P 
function is given by a convergent power series (as in all examples we examine here), 
then the resulting series for x{Q'^/fi'^,XQ) has a finite radius of convergence at any 
fixed \Q'^/fi^\, and defines the unique correct analytic continuation to complex Q^. 
Given the finite convergence radius, there are no singularities for a fixed IQ^/a'-^I if 
Xq is small enough. Due to asymptotic freedom, small enough xq at fixed IQ^/a*^] 
corresponds to large \Q'^\ and therefore x{Q^) has no singularities for large enough 
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Considering the 2-loop APT spectral function, PR{t) — —{l/TT)lm{x{—s — ie)} 
with x{Q'^) given by eq. (11), we show in fig. 4 the vahics of W along the time-like 
axis, below the cut, for a 2-loop f3 function with /5o = 1 and several different values 
of c. We identify two categories of lines in fig. 4: 

(a) for c < —f3o (c < — 1 in the figure) W fiows to zero in the infrared. This leads 
to a non-trivial perturbative infrared fixed-point: x{Q'^) — > xpp = — 1/c. 

(b) for c > — /?o (where c can be either positive or negative) \W\ — > oo in the 
infrared, and thus xiQ"^) — > 0. 

When c < — /^o, the singularity structure (see fig. 3) is consistent with causality: 
the cut is only on the time-like axis. This guarantees that x{Q'^) obeys the dispersion 
relation in eqs. (5) and (6). It follows that in this case the APT coupling coincides 
with the perturbative coupling: xpji^iQ"^) = x{Q'^). On the other hand, when 
c > —Po, the singularity structure is always inconsistent with causality (see fig. 3). 
Therefore in this case xjsj3j^{Q'^) differs from the perturbative coupling x{Q^). 

Having identified how the relative size of c vs. affects the main features of the 
APT coupling, we are now ready to consider the actual 2-loop f3 function in QCD 
and examine the dependence of the analytically continued coupling x{—s) on Nf, 
the number of light quark flavors. In fig. 5 we plot values of W along the time-like 
axis, below the cut, for Nf = 0, 1, 2, . . . , 16. As already mentioned, the branch Wn 
in which the solution corresponding to the time- like axis resides, depends on c/j3o- 
The table below lists the values of Nf, and the corresponding branch labels n. 



Nf 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


n 


1 


1 


1 


1 


1 


1 


1 


2 


25 


-1 
























For A^/ > 10 we have c < — /3o, so the APT couphng coincides with the perturbative 
one. On the other hand, for A^j < 9 we have c > — /^o, and then we expect the APT 
solution to be different. For Nf = 9, —(3o < c < and therefore a pair of complex 
singularities appears, while for Nf < 8 there is a single singularity of the space-like 
axis. The Nf = 8 line is not shown in fig. 5, since the large value of Po/c in this 
case causes the corresponding Im(iy) to fall within the W25 branch, way outside the 
vertical range of the plot. It is important to note that for the physically relevant 
case, Nf < 3, the exact number of fiavours has only a small effect on W. 

Next, we use (11) to calculate x{—s) below the cut for two representative exam- 
ples: Nf = 3 and Nf = 10. The corresponding real and imaginary parts of x{—s) 
are shown in fig. 6 as function oi t — ln(s/A^). The ultraviolet limit of the real 
part is the same in both cases: Re{x(— s)} — > 0. The difference is in the infrared: 
for Nf = 10, Re{x(— s)} xpp = —1/c, in accordance with (a) above, while for 
Nf = 3, Re{a;(— s)} 0, in accordance with (b) above. The imaginary part of 
x{—s) vanishes in both the infrared and in the ultraviolet limit, so that the APT 
coupling given by (5) and (6) is well defined. 

Im{x(— s)} below the cut is always positive-definite and thus the APT coupling 
is a monotonically decreasing function of Q^, attaining its maximum value a^APTlO) 
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at = 0. The infrared limit of the APT couphng can be obtained by integrating 
Im{x(— s)} over all scales: 

^APt(O) = - / lm{x{-s)} dt (27) 

This integral can be performed analytically by changing the integration variable to 
W. Using (10) and (11) we write dt = (-c//3o)(l + W)/WdW and obtain: 



11 r /-uv 1 

.,P,(0)^--In,{/^ -} (28) 



where the integration contour in the W plane is along the line that corresponds to 
the time-like axis (examples are provided by the continuous lines in figs. 4 and 5). 

For c > and for — /3o < c < the time-like axis is always mapped into a non- 
principle branch (n ^ 0) and then the integration contour stretches from Re{l^} — 
— oo to Re{VF} — > +00, while lm{14^} is bounded. The integral can then be 
performed by closing the integration contour with a semi-circle at infinity. Since 
there are no poles inside the closed contour and the integration along the semi-circle 
yields ztt, one obtains xp^^i^) = I/Pq. 

For c < —Po the relevant branch is n = and then the contour starts at 1^ = in 
the infrared and reaches Re{VF} — > -\-oo in the ultraviolet. In this case we evaluate 
the integral directly and then use the asymptotic behavior of Wo{z) at small and 
large \z\. The result is a^ApxlO) = — 1/c. A posteriori, this should not come as a 
surprise, since we already know that whenever the time-like coupling is within the 
Wq branch, (c < — /3o), the APT solution coincides with the perturbative couphng 
which fiows to the infrared fixed-point at —1/c. 

We note that the conclusion of refs. [4, 5] that the infrared limit of the 2-loop 
APT coupling is always l/f3o is correct only for c > —(3o- 

Our results for 2-loop a;APT(0) summarized in fig. 7, where the upper plot 
(continuous line) shows the 2-loop APT infrared hmit for a /3 function with Po — 1 
and a span of c values, and the lower plot (continuous line and crosses) shows the 
2-loop APT infrared limit for QCD with a varying number of flavours. For Nf < 9 
the APT infrared limit is l//9o, while for Nf > 10 it is —1/c. 

In a similar manner we now calculate the infrared APT limit for the PA-improved 
three-loop coupling defined from the imaginary part of x{—s) in eq. (22). Also here 
we perform the integral analytically by changing the integration variable to W, 

1 1 r /-uv w + 1 dW] 
,^p,(0) = --to|/^ _____|. (29) 

We find the same two scenarios: when the time-like axis corresponds to a non- 
principal branch, we close the contour by a half-circle at infinity. Otherwise, for 
n = 0, we evaluate the integral directly and use asymptotics of Wo{z). 

Depending on the values of /9o, c and C2, there are four different possibilities, as 
summarized below: 
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3^apt(0) 



c> 
c < i 



( c>-(3o 
c < -Po 



C2 < 
C2 > 



l//5o 
l//9o 

(-1/C)/(1 - C2/C^) 

l//3o-(l/c+l//?o)/(l-C2/c2) 



(30) 



For c > — /5o we again obtain a completely universal infrared limit: a;APT(0) 
does not depend at all on c and C2. The universality breaks down for c < — /?o, 
just like in the two-loop case. When the singularity structure of the perturbative 
coupling is consistent with causahty, we expect that it coincides with the APT 
couphng, and then xjsjyj-{0) should be just equal to the solution of f3{x) = 0, which 
is = 0) = (— l/c)/(l — 02/0^). However, contrary to the two-loop case, at 

three-loops, given c < —/So, it is still possible that the singularity structure will be 
inconsistent with causality: if C2 > the perturbative coupling has an extra pole 
on the space-like axis (see discussion below eq. (23), case (b)). Indeed, in this case 
the APT infrared limit is different from the perturbative one. 

Our results for the infrared limit of the APT coupling at 3-loops are shown in 
fig. 7, on top of the 2-loop results. The plot demonstrates that as long as c > — /3o the 
limit a^APT(O) = l/Z^o is universal, i.e. it does not depend on c and C2. Nevertheless, 
when c < —Pq the values of c and C2 are important for the APT infrared limit, as 
indicated by (30). 

In the upper plot, one notes that for a 3-loop case with C2 > and c < 0, 
^APt(O) can become large, and even diverge, when C2 ~ c^. Then a; apt ceases to 
be a good expansion parameter. Another interesting feature of the 3-loop solution 
is the jump that may occur in the value of a;APT(0) at c = — /3o when c is varied. 
Such a jump indeed occurs if C2 < (dashed line in the upper plot). 

In the lower plot, one observes that the APT coupling based on the PA-improved 
3-loop P function in MS shows a sharp transition at Nf ~ 10. For Nf > 10 the 
coupling has a perturbative infrared fixed point at a rather small value (a;(0) ~ 0.1), 
suggesting a reliable purely perturbative phase. On the other hand, for Nf < 9 
the APT coupling has a non-perturbative infrared fixed point, with the "universal" 
value of l/Po- It is natural to ask whether the transition observed at Nf ~ 10 
is an indication of the phase transition expected in QCD as Nf is increased (see 
[3] and refs. therein), or is it an artifact of considering a truncated /3 function in 
a special renormalization scheme. While we cannot fully answer this question, it 
is worthwhile noting that the transition point Nf ~ 10 docs not depend on the 
renormalization scheme, since it refiects the scheme-independent condition c = — /9o- 
Moreover, the essential qualitative features of the MS results for Nf > 10 are likely 
shared by most acceptable renormalization schemes: one observes that in physical 
schemes the condition C2 < is always obeyed for Nf > 10 (see fig. 1 in [3]). 
Therefore, the perturbative fixed-point of the PA-improved 3-loops /3 function at 
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( — 1/c) / ( 1 — C2/c^) is positive and the perturbative coupling has a causal analyticity 
structure. The last case considered in (30), i.e. c < — /3o and C2 > c^, is probably 
never realized in QCD. 

6 Summary and Conclusions 

Using the Lambert W function we achieved a thorough understanding of the struc- 
ture of the solutions of the 2-loop and (Pade improved) 3-loop RG equation for the 
coupling constant in the complex plane. The main result is that the running 
perturbative coupling in QCD can have a causal analyticity structure, with a non- 
trivial infrared fixed point, provided /9o is small enough and that i3\ is negative, 
i.e. in the framework of the Banks-Zaks expansion. This suggests that a consis- 
tent, purely perturbative definition of QCD should be possible when the number of 
fiavors is large enough (barring possible complications due to renormalons and to 
the asymptotic nature of the (5 function series and of the Banks-Zaks series, which 
have not been addressed here). On the other hand, for larger /Sq or for positive 
unphysical singularities are present in the infrared region, which are to be removed 
by non-perturbative power-like effects. 

At the 2-loop level, the causality requirement translates into the condition 
/9i < — /3q, i.e. xpp = — 1/c < l//3o- This condition is also relevant for the PA- 
improved 3-loop case, for most of the acceptable renormalization schemes. Taking 
seriously the 2-loop and 3-loop results, we are lead to conclude that a purely per- 
turbative phase is reahzed in QCD with Nf > 10. This result is in reasonable 
agreement with other approaches (sec [3] and refs. therein). 

We have investigated the APT approach, a simple mathematical way (which 
admittedly lacks a physical basis, see the discussion below) to implement the nec- 
essary power-like effects, starting from perturbation theory. In this approach, a 
causal coupling is reconstructed via a dispersion relation from the time-like discon- 
tinuity of the perturbative coupling. We have shown that in some cases the value 
of the APT infrared fixed point may depend on higher-order coefficients of the per- 
turbative f3 function, and is not always given by the one-loop value l//3o, which is 
therefore not "universal". Departure from the "universal" value appears not only 
in the above mentioned cases where the perturbative coupling is by itself already 
causal and the APT fixed point coincides with the perturbative one, but also when 
unphysical space-like singularities are present and a (positive) perturbative infrared 
fixed point does not exist - case (b) below eq. (23). Nevertheless, the latter case 
might be viewed as academic, since in QCD it is not realized for any number of 
flavors. 

It is therefore natural to wonder what general conditions are required to re- 
cover the "universal" APT value, which clearly has a special, remarkable status (in 
particular, all the "non-perturbative" APT fixed point curves end up on the "uni- 
versal" l//3o line in Fig. (7)). Sufficient conditions for the "universal" APT value 
were given in [7], assuming a Landau singularity is present on the space-like axis: 
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essentially, the perturbative coupling was required to approach the trivial infrared 
fixed point when the scale is decreased to zero on the space-like axis below the 
Landau singularity**. Thus, in order to avoid universality in this case, there must 
exist a non-trivial, but perturbative, infrared fixed point (at unphysical negative or 
complex values of the coupling), whose domain of attraction includes the trajectory 
going through the Landau singularity. Indeed, a (negative) perturbative infrared 
fixed point is present in the 3-loop example for c < and C2 > c^, explaining why 
the last "non-perturbative" APT fixed point in eq. (30) differs from the "universal" 
value. Actually, universality may be obtained under more general conditions then 
those stated in [7], e.g. in the case — /?o < c < in eq. (30), where the APT value 
is still the universal one, despite the presence of a non-trivial perturbative infrared 
fixed point (which always exists for c < 0), and regardless of its sign (if the latter is 
negative, then a Landau pole is still present on the space-like axis). 

In all the examples presented here, namely the 2-loop and PA-improved 3-loop 
(3 functions, universality is obtained whenever the perturbative coupling approaches 
the trivial infrared fixed point as — > 0" on the time-like axis. This has a simple 
mathematical explanation in terms of the Lambert W function, as in these cases, the 
contour in the W plane stretches from Re{VF} — > —00 to Rc{VF} — > +00, while 
Im{iy} is bounded. The integral can then be performed by closing the integration 
contour with a semi-circle at infinity. Since there are no poles inside the closed 
contour and the integration along the semi-circle yields iir, one obtains a;APT(0) — 
1/Po- As a result, the details of the contour, which do depend on the coefficients of 
the f3 function, are insignificant. 

It would be interesting to know whether, for a generic higher-order (3 function, 
a universal 1/Po limit of the APT coupling is obtained whenever the perturbative 
coupling approaches the trivial infrared fixed point on the time-like axis, giving a 
sufficient condition for universality alternative to that of [7]. Unfortunately, because 
of the complicated phase structure in the complex plane when complex Landau 
singularities are present (see the discussion before eq. (17)), it is not clear if such 
a condition could turn out to be also necessary. Note that in the present examples, 
for — /3o < c < there is a 2-phase structure, such that the space-hke couphng 
approaches the non-trivial perturbative infrared fixed point, while the time-like cou- 
pling approaches the trivial one. On the other hand, for c < —f3o < there is only 
one phase and both couplings approach the non-trivial fixed point. But, in principle, 
there could be more complicated examples where a third scenario is realised, namely 
the space-like coupling approaches the trivial fixed point (thus insuring universality 
according to [7]), while the time-like coupling approaches a non-trivial one. 

It is also interesting to note that the PA-improved 3-loop coupling probably 
offers the simplest example (in case c < and C2 > c^) where the standard relation 
between renormalons and Landau pole (sec previous footnote) might not hold, since 
the Landau pole trajectory is not in the domain of attraction of the trivial infrared 
fixed point in this example. Developing this point is beyond the scope of the present 

**This condition also allows to relate [21] the ambiguity arising from integrating over the Landau 
singularity to the one due to renormalons, a fact exploited in the proof. 
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paper. 

One might question the physical relevance of the APT coupling. It is unlikely, 
given its perturbative origin, that it exhausts all non-perturbative effects in the 
running coupling, but it may play a useful phenomenological role (yet to be clarified) 
in a more general framework [7, 22] . 

Acknowledgments 

This research was supported in part by the Israel Science Foundation administered 
by the Israel Academy of Sciences and Humanities, by a Grant from the G.I.F., the 

German-Israeli Foundation for Scientific Research and Development, by the Charles 
Clore doctoral fellowship, and by the EC program 'Training and Mobility of Re- 
searchers', Network 'QCD and Particle Structure', contract ERBFMRXCT980194. 

References 

[1] D. Gross and F. Wilczek, Phys. Rev. Lett. 30 (1973) 1343; H.D. Politzer Phys. 
Rev. Lett. 30 (1973) 1346. 

[2] V. I. Zakharov, Renormalons as a bridge between perturbative and nonpertur- 
bative physics, hep-ph/9802416. 

[3] E. Gardi, M. Karliner, Relations between observables and the infrared fixed-point 
in QCD, to be pubhshed in Nucl. Phys. B, hep-ph/9802218. 

[4] I.L. Solovtsov and D.V. Shirkov, Phys. Rev. Lett. 79 (1997) 1209-1212, An- 
alytic model for the QCD running coupling with universal as{0) value hep- 
ph/9704333. 

[5] D.V. Shirkov, On the analytic "Causal" model for the QCD running coupling, 
hep-ph/9708480; K. A. Milton, O. P. Solovtsova, OKHEP-97-06, Analytic per- 
turbation theory: a new approach to the analytic continuation of the strong cou- 
pling constant alphas into the time-like region, hep-ph/9710316; I.L. Solovtsov 
and D.V. Shirkov Analytic approach to perturbative QCD and renormalization 
scheme dependence, hep-ph/9711251; K. A. Milton, I.L. Solovtsov and O. P. 
Solovtsova, Phys. Lett. B415 (1997) 104. 

[6] Yu.L. Dokshitzer, G. Marchesini and B.R. Webber, Nucl. Phys. B469 (1996) 
93. 

[7] G. Grunberg, Power corrections and Landau singularity, hep-ph/9705290. 

[8] R. Akhoury and V.I. Zakharov, UM-TH/97-11, hep-ph/9705318; UM-TH/97- 
19, hep-ph/9710257; UM-TH/97-18, hep-ph/9710487; G. Grunberg, CERN- 
TH/97-340, hep-ph/9711481. 



17 



[9] W. E . Caswell, Phys. Rev. Lett. 33, (1974) 244; D. R. T. Jones, Nucl. Phys. 
B75, (1974) 531. 

[10] T. Banks and A. Zaks, Nucl. Phys. B196 (1982) 189. See also G. Grunberg, 
Phys. Rev. D46, (1992) 2228. 

[11] R.M. Corless, G.H. Gonnet, D.E.G. Hare, D.J. Jeffrey and D.E. Knuth, "On the 
Lambert W function". Advances in Computational Mathematics, 5 (1996) 329, 
available from http : / /pineapple . apmaths . uwo . ca/~ rmc/papers/LambertW/ 

[12] C. J. Maxwell, Phys. Lett. B409 (1997) 450-454. 

[13] G. Cvctic, Phys. Rev. D57 (1998) 3209-3213; G. Cvctic, Nucl. Phys. B517 
(1998) 506-520; see also G. Cvetic and R. K5gerler, MPI-PliT 98-13, BI- TP 
98/2, hcp-pli/9802248. 

[14] E. Gardi, Phys. Rev. D56 (1997) 68; S. J. Brodsky, J. Ellis, E. Gardi, M. 
Karliner, M.A. Samuel. Phys. Rev. D56 (1997) 6980-6992. 

[15] G. Grunberg, Phys. Rev. D29 (1984) 2315. 

[16] P.M. Stevenson, Phys. Rev. D23 (1981) 2916. 

[17] N.N. Khuri and H.C. Ren, Annals Of Physics 189 (1989) 142-154. 

[18] G. Grunberg, Phys. Lett. B304 (1993) 183. 

[19] N.G. Uraltsev, private communication with G. Grunberg. 

[20] M.A. Samuel, J. Ellis and M. Karhner, Phys. Rev. Lett. 74 (1995) 4380; Phys. 
Lett. B400 (1997) 176; J. Ellis, E. Gardi, M. Karliner and M.A. Samuel, Phys. 
Lett. B366 (1996) 268; Phys. Rev. D54 (1996) 6986; I. Jack, D.R.T. Jones and 
M.A. Samuel, hep-ph/9706249; J. Elhs, I. Jack, D.R.T. Jones, M. Karliner and 
M.A. Samuel, Phys. Rev. D57 (1998) 2665. 

[21] G. Grunberg Phys. Lett. B372 (1996) 121. 

[22] A.I. Alekseev, QCD Running Coupling: Freezing Versus Enhancement in the 
Infrared Region, hep-ph/9802372; A.I. Alekseev, B.A. Arbuzov, Analyticity and 
Minimality of Non-perturbative Contributions in Perturbative region for as, 
hep-ph/9704228. 



18 




-10' — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ' 

-20 -10 10 20 

Re(W) 



Figure 1: The complex W plane is divided by dashed lines to branches of the 
Lambert W function: Wq, W±i, W±2- The dashed line that forms the boundary 
between Wi and represents also the values of W for the space-like coupling at 
2- loops for c > (above the Landau singularity). The continuous lines correspond 
to the values of W for a set of fixed phase rays in the lower half-plane for the 
example of a 2- loop P function with /9o = 1 and c = 2/1. Large corresponds 
to the left side of the lines, while the infrared limit corresponds the right side. The 
phases in the plane arc (f) = —gn. where g = 0.1, 0.2, 0.3, 0.9. The dot- 
dashed line {g = 1) corresponds to the time-like axis below the cut. The use of W2 
for (f) < —01 = —2{c/Po)7r ~ — O.STItt guarantees the continuity of W. 
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Figure 2: The complex W plane is divided by dashed lines to branches of the 
Lambert W function: Wq and W±i. The g = line corresponds to the values of 
W for the space- like coupling at 2-loops for c < 0. The rest of the continuous lines 
correspond to the value of W for a set of fixed phase rays in the lower Q'^ half-plane 
for the example of a 2-loop /3 function with Pq — 1 and c = —2/3. Large \Q'^\ 
corresponds to the right side of the lines. The phases in the plane are = —qti, 
where g = 0.1, 0.2, 0.3, 0.9. The dot-dashed line {g = 1) corresponds to the 
time-like axis below the cut. The use of W-i for (f) < — 0o = ~|c//3o|7r ~ — 0.6677r 
guarantees continuity of W for > Ql (see eq. (16)). 
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Figure 3: The three possibihtics for the singularity structure of the 2-loop couphng. 
Branch points are represented by diamonds and cuts by bold lines. Only for c < —f3o 
(upper plot) perturbative freezing leads to a singularity structure that is consistent 
with causality. For — /3o < c < there is a pair of complex conjugate branch points 
and for c > there is one space-like branch point that violate causality. Note that 
with the present analytical continuation, the non-causal singularities are confined 
to the infrared region, where they can be removed from physical quantities by non- 
perturbative effects. 
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Figure 4: The various continuous lines show the value of W along the time-like axis 
(below the cut), corresponding to the 2- loop coupling for = 1 and different values 
of c. For the c > large IQ^I corresponds to the left side of the lines, while for 
c < large \Q^\ corresponds to the right side. The dashed lines show the division 
of branches of the Lambert W function: Wq, W±i and W±2- 
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Figure 5: The various continuous lines show the value of W along the time-like axis 
(below the cut), corresponding to the QCD 2- loop coupling for different value of 
iVj = 0,l,2,...,16. The dashed lines show the division of branches of the Lambert 
W function: Wq and W±i. 
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Figure 6: The real (dashed hne) and imaginary (sohd hne) parts of the analytically 
continued coupling on the time-like axis, x{—s), for the case of a 2-loop (3 function 
in QCD with A^/ = 3 % = 9/4 and c = 16/9) in the upper box, and Nf = 10 
{Po — 13/12 and c = —37/26) in the lower box. For Nf — 10 perturbativc freezing 
leads to a well defined infrared limit for the perturbative coupling: x^^ — — 1/c, 
while for Nf — 3 the "perturbative infrared limit" is zero. 
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Figure 7: The upper box shows the infrared hmit of the APT solution for a hy- 
pothetical PA-improved 3- loop /3 function with /Sq — 1, a span of values for c, 
—12 < c < 12 and various values for C2- Continuous line: C2 = (the 2-loop case), 
dashed line: C2 = —1.4, dot-dashed line: C2 = 20, dotted line: C2 = 200. The lower 
box shows a;A.PT(0) QCD with < Nf < 16. Squares and crosses correspond 
to a 2-loop and a PA-improved 3- loop MS (3 function, respectively. For Nf < 9, 
^APt(O) is l//9o in both cases. 
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